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ABSTRACT 


An  investigation  of  the  ignition  delay  of  hydrogen- air  mixtures  near  the  second  explosion 
limit  at  pressures  of  15  and  30  psia  was  made  in  a  shock  tube.  The  shock  tube  calculations 
are  discussed  along  with  the  problems  encountered  during  the  investigation. 

Attempts  to  correlate  the  data  were  satisfactorily  accomplished  when  the  H09  recombi- 

-0  4  * 

nation  .rate  was  assumed  to  be  proportional  to  T  V. 
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SECTION  I 

INTRODUCTION 


The  mechanism  of  the  hydrogen-oxygen  reaction  has  been  studied  for  many  decades.  In 
spite  of  the  simple  overall  reaction,  the  actual  mechanism  has  been  subject  to  controversy 
for  a  number  of  years.  If  a  reaction  scheme  is  known,  along  with  the  rate  constants  for  these 
reactions,  the  composition  and  temperature  of  the  hydrogen-oxygen  mixture  can  be  computed 
as  the  mixture  reacts.  Ignition  delay  can  be  defined  in  a  number  of  ways,  but  the  definition 
must  be  such  that  the  delay  is  experimentally  detectable  in  order  to  check  the  calculations. 
Two  of  the  most  prevalent  definitions  of  ignition  delay  are:  (1)  that  point  in  time  when  there 
is  a  sudden  increase  in  the  pressure  of  the  mixture  and  (2)  that  point  in  time  when  there  is 
a  sudden  increase  in  the  concentration  of  OH. 

Various  methods  of  experimentally  determining  ignition  delays  have  been  used.  The  principal 
ones  Include  well-stirred  reactors,  flowing  systems,  and  shock  tubes. 

The  well-stirred  reactor  is  a  closed  vessel  containing  heated  air  into  which  hot  hydrogen  is 
introduced  and  rapidly  mixed  with  the  air.  The  pressure  of  the  vessel  is  monitored  for  a  sudden 
rise,  wnich  signifies  the  end  of  the  ignition  delay  period.  The  two  disadvantages  of  this  method 
are  that  the  mixing  tends  to  obscure  the  results  when  ignition  delays  on  the  order  of  a  few 
hundred  microseconds  are  being  considered,  and  the  results  tend  to  depend  on  the  vessel  size 
and  wall  coating. 

In  the  second  method,  hydrogen  is  introduced  into  a  flowing  system  of  subsonic  air  at  a 
temperature  slightly  below  the  ignition  limit.  A  diverging  section  slows  the  air  down  and  raises 
the  temperature  of  the  mixture  above  the  ignition  limit.  The  distance  from  the  diverging  section 
to  where  combustion  occurs,  or  the  combustion  zone,  is  used  to  calculate  the  ignition  delay. 
With  this  method,  however,  it  is  difficult  to  obtain  a  uniform  mixture  and  we  must  assume 
that  no  reaction  has  occurred  before  the  mixture  enters  the  diverging  section. 

The  third  method  has  become  popular  with  the  advent  of  the  shock  tube.  The  hydrogen-air 
mixture  cbntained  in  the  driven  section  of  the  shock  tube  is  heated  rapidly  to  the  correct  tem¬ 
perature  and  pressure  by  a  normal  shock  wave.  The  time  interval  elapsed  from  the  instant  of 
the  passage  of  the  shock  wave  to  the  moment  when  the  OH  radiation  appears  is  measured. 
When  measurements  are  made  behind  the  incident  shock,  a  diluent,  usually  argon,  is  added  to 
the  hydrogen-oxygen  mixture  to  prevent  the  heat  released  by  the  reaction  from  forming  a 
detonation  wave  behind  the  normal  shock  wave.  Measurements  behind  the  reflected  shock 
wave  are  also  affected  by  the  tendency  to  form  detonation  waves,  and  measurements  cannot 
be  made  if  the  reflected  shock  wave  is  not  past  the  measuring  station  before  the  heat  is  re¬ 
leased.  In  this  study,  ignition  delays  on  the  order  of  80-100  microseconds  were  measured 
before  the  formation  of  detonation  waves  became  a  problem. 

Shock  tube  studies  of  the  ignition  delay  in  mixtures  of  hydrogen-oxygen  have  given  support 
to  a  postulated  reaction  scheme  and  indicate  that  the  Initiation  reaction  is 

H2  +  0?  —  H2  0* -  2 OH 


and  not 


H2  — —  2  H 

Although  the  reaction  scheme  does  predict  correct  values  of  the  ignition  delay  at  temper¬ 
atures  above  1100-1200°K,  it  does  not  properly  account  for  the  occurrence  of  the  Becond 
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explosion  limit  of  the  hydrogen-oxygen  reaction.  To  account  for  this  limit,  it  was  postulated 
that  the  reaction 


H  +  02  -4-  M  — —  H02  +  M 

leads  to  a  termination  of  the  chain.  HOg  has  never  been  observed,  however,  so  this  assump¬ 
tion  is  not  entirely  satisfactory,  although  a  paramagnetic  resonance  spectra  for  the  H02 

molecule  was  observed  recently  (Reference  1)  and  a  rate  constant  for  the  above  reaction  was 
established. 

With  this  additional  reaction,  analytical  expressions  for  the  ignition  delays  of  hydrogen- 
oxygen  mixtures  have  been  derived  (Reference  2).  These  results  explain  the  long  ignition 
delays  as  the  second  explosion  limit  Is  approached.  Data  in  this  regime  are  scarce;  only 
recently  a  comprehensive  study  was  made  (Reference  3),  but  these  observations  do  not  seem 
to  follow  any  general  pattern  as  to  pressure  dependence  and  they  differ  from  the  values  ob¬ 
tained  by  the  analysis  of  Reference  2.  Therefore,  there  Is  some  doubt  as  to  the  accuracy  of 
these  data. 

On  the  basis  of  this  information,  it  was  decided  to  Btudy  the  Ignition  delay  of  hydrogen- 
air  mixtures  under  carefully  controlled  conditions.  Near  the  second  explosion  limit,  the 
ignition  delay  time  changes  drastically  with  small  changes  in  temperature;  therefore,  the 
absolute  values  of  ignition  delay  cannot  be  determined  accurately.  The  general  trend,  however, 
is  expected  to  provide  a  check  on  the  analysis  of  Reference  2. 

The  reflected  shock  technique  was  used  in  this  investigation  so  that  we  would  not  have  to 
dilute  the  hydrogen-air  sample  to  prevent  the  formation  of  detonation  waves  behind  the  in¬ 
cident  shock  when  the  ignition  delay  was  short.  In  addition,  the  test  time  behind  the  incident 
shock  was  not  sufficient  for  the  case  when  the  ignition  delays  were  long.  The  tailored  inter¬ 
face  mode  of  operation  was  also  used  since  test  times  as  long  as  10,000  microseconds  might 
be  required.  In  the  tailored  interface  mode,  the  driver  gas  composition  is  adjusted  such  that 
the  reflected  shock  wave  will  pass  through  the  contact  surface  without  propagating  any  dis¬ 
turbance  into  the  test  gas. 
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section  n 

SHOCK  TUBE  CALCULATIONS 


Since  accurate  temperature  measurements  behind  the  reflected  shock  are  not  possible  at 
present,  one  must  usually  calculate  the  temperature  from  the  measured  shock  wave  velocity 
and  initial  conditions.  From  Reference  4,  the  general  basic  equations  for  a  normal  shock  wave 
are: 


Continuity: 


p  W  --P  [w  ±  Cu,  -u2)] 


Momentum: 


Energy: 


Pi  2 

p  4_  wz  =  p 

i  q  2 


h.  4-  — 


+  -A-  [  W*  (.,  -u, l] 


where  the  (+)  sign  refers  to  a  right-traveling  wave  and  the  (-)  sign  refers  to  a  left-traveling 
wave. 

For  the  general  case  of  a  reacting  gas  or  a  gas  with  variable  specific  heats,  these  equations 
must  be  solved  by  an  iteration  procedure.  The  conditions  behind  the  incident  shock  are  obtained 
by  combining  Equations  1,  2,  and  3. 

Il  |  +  <"i-M  Jj_  (4, 

p.  11+  p,  /  >  p, 

One  first  specifies  a  value  for  and  then  assumes  a  value  for  P2  (obtained  from  an  ideal 
gas  formula).  From  this  value  of  P 2  and  T^,  h2  and  p^  can  be  calculated  and  a  new  value  for 
P0  calculated  from  Equation  4.  This  calculated  value  of  P2  then  becomes  the  new  assumed 
value  of  P2  and  the  process  is  repeated  until  the  assumed  and  calculated  values  of  P2  are  the 
same.  The  shock  wave  velocity  is  then  computed  from 

« .  /Tr~rp,Tp,~-rr7  , 

V  P ,  (I -  p  /  p  ) 

I  2 

The  particle  velocity  behind  the  incident  shock  is 


(  P,  i  Pg  A 

f.  P'  '■  „ 

~r  i  p  "  1 ) 
p.  \  n  > 

V'TD  9 

where  is  normally  zero. 
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For  the  reflected  shock,  the  process  becomes  a  little  more  complicated  since  there  are 
two  equations  upon  which  to  iterate.  These  equations  are 


+  K.hg  +  -L  -fLJ-LftJ.] 


(pg/p.-iHI-fl  /pz)  - 

f  1  _  l  (Pg/P.  >  Ipt'pj)  J  1 
^  1 "  P9  /  Pa  - 1  * 


One  first  assumes  a  value  for  Pg  from  ideal  gas  formula  and  calculates  a  value  for  Tg  from 
Equation  8  by  assuming-^.  equals  With  this  value  of  Tg  and  Pg,  ^  can  be  calculated  and 
a  new  value  of  Tg  computed.  This  is  repeated  until  two  successive  Tg's  are  the  same.  A  value 
for  Pg  is  then  computed  from  Equation  7  and  the  entire  process  is  repeated  until  two  succes¬ 
sive  values  of  Pg  are  the  same.  The  reflected  shock  velocity  is  then  found  from 

“u2 

urK  =u2-W=- - s- - -  (9) 


[l-Pz'p9) 


For  the  simplest  case  when  operating  in  the  reflected  shock  mode,  the  test  time  is  normally 
interrupted  by  an  expansion  wave  or  shock  generated  when  the  reflected  shock  passes  through 
the  discontinuity  (contact  surface)  between  the  driver  and  driven  gases.  This  test  time  can  be 
extended  by  operating  the  shock  tube  in  the  tailored  mode,  which  can  be  accomplished  by  one 
of  two  methods.  The  first  method  is  to  heat  the  driver  gas  to  the  proper  temperature,  which 
limits  operation  to  shock  speeds  greater  than  the  tailoring  shock  Mach  number  of  the  unheated 
driver  gas.  The  second  method  limits  operation  to  shock  speeds  below  that  for  the  pure  driver 
gas:  this  is  accomplished  by  introducing  a  foreign  gas  into  the  driver  gas,  thus  altering  the 
molecular  weight  and  ratio  of  the  specific  heats  of  the  driver  gas. 

To  calculate  the  conditions  required  for  tailored  operation,  one  first  writes  the  equations 
of  the  difference  in  particle  velocity  across  the  reflected  shock  in  the  driver  and  driven  gases: 


u3  "  U2 


P2  (P3/P2  -I  )  (I  -pz  /  p9)  /Pz 


UT  ~U3  =  V  P3  <P7  /P,  *  I  5  (  I  ~  p3  /  p7  )  t  p3 
Now  Pg  =  Pg,  Ug  *  Ug,  a nd  Ug  -  0.  The  conditions  required  for  tailoring  are  that  P7  =  Pg  and 
u?  =  0. 

One  then  arrives  at  the  conditions  required  for  tailoring  as- 

„01 

which  can  be  rewritten  as 


;p9/p2-  (»T2/fre)(Tg/T2)j _ 

>,/P2)  (»?,/»,)-  (W2  /Wy)(T/T,) 
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For  the  case  of  an  ideal  gas  in  the  driver  section,  Equation  11  can  be  written  as 

li  -  [Pg/Pj-  <t,/t2)  ]  [i+a4  <pg/pt)] 

t2  '  (m2  /wA)  (Ps  /p* )  (p8/p2'I)  (a4*i )  11 

Since  the  reduction  of  conditions  in  the  driver  gas  is  the  result  of  a  series  of  infinitesimal 
isentroplc  expansion  waves,  Tg  may  also  be  written  as 

f  1Z 

t*  -  t< 

To  calculate  the  driver  gas  composition  required  for  tailoring  at  shock  speeds  below  the 
tailoring  shock  speed  for  the  pure  driver  gas,  one  assumes  a  composition  for  the  driver  gas 
and  calculates  from  Equations  12  and  13.  New  compositions  are  assumed  until  from 

both  equations  agree. 


The  ideal  driver  pressure  can  be  calculated  from 


vp, 


2a 


(14) 


For  a  given  shock  speed.  Equation  14  yields  only  an  approximation  of  the  pressure  ratio  re¬ 
quired  to  produce  a  shock  of  given  strength  since  the  boundary  layer  growth  behind  the  shock 
wave  attenuates  the  shock  as  it  proceeds  away  from  the  diaphragm  station. 


Equations  4  through  14  were  programmed  for  the  IBM  7094  to  calculate  all  the  conditions 
in  the  Hg-air  test  gas  for  shock  Mach  numbers  between  2.2  and  3.2;  the  calculations  allowed 

for  variable  specific  heats  but  no  chemical  reactions.  Another  program  was  also  written  to 
calculate  the  conditions  behind  the  shock  wave  for  air,  allowing  the  gas  to  bo  in  chemical  and 
thermodynamic  equilibrium  but  assuming  no  ionization,  for  use  during  initial  firing  and  cali¬ 
bration  of  the  shock  tube. 


Figures  1  and  2  show  the  reflected  shock  temperature  and  pressure  and  were  used  to  de¬ 
termine  the  test  conditions  from  the  measured  shock  speed.  For  long  delay  times  (greater 
than  2  milliseconds),  tailored  operation  is  required.  Figure  3  was  used  to  determine  the 
amount  of  nitrogen  that  had  to  be  added  to  the  helium  driver  for  tailored  operation. 

Figure  4  shows  the  required  pressure  ratio  across  the  diaphragm  to  produce  a  shock  of 
given  strength.  Experimental  pressure  ratios  are  given  to  indicate  the  degree  of  shock  wave 
attenuation.  Figure  5  shows  the  pressure  behind  the  incident  shock.  The  lowest  allowable 
value  of  reflected  shock  pressure  was  used  in  these  tests  to  produce  an  incident  shock  pres¬ 
sure  that  was  great  enough  to  trigger  the  electronic  instrumentation. 

Figure  6  indicates  the  temperature  behind  the  incident  shock,  which  determines  the  upper 
limit  of  shock  speed  that  can  be  investigated,  since  Tg  must  be  low  enough  that  no  reactions 

of  the  Ilg-air  mixture  occur  behind  the  incident  shock. 
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Figure  2.  Reflected  Shock  Pressure  Versus  Shock  Mach  Number 
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SECTION  III 

IGNITION  DELAY  COMPUTATIONS 

In  Reference  2  the  following  reactions  were  assumed  to  be  of  importance  during  the  induc¬ 
tion  period: 

OH  +  Hg  -^HaO+H 
k  a 

H  +  Og  — •  OH  +  0 

0  +  H  a  OH  ♦  H 

k« 

H  +Of  +  M—i*HO*+M 
HO*  +  H,-^  H20t  +  H 

where 


kj  *  6.3  *  10 10  e  -S900/RT  .0/mole-sec 

k2  =  4.0  x  IOn  #-^000/RT  .0/mole- sec 

k3=  1.2  x  I0H  •  "W50/RT  0 /mol* -sec 

k-s  3.27  x  I015  (  xH  +  0.35  xft  +  0.43  xu  +  0.2  x. 

•  "a  Oj  Nj  Ar 

2  2 

+  14.3  xM^0  + - )  T”*'92  t  /mole-sec 

k„=  5.4  x  10*  e'24,000/RT  i/mol# -sec 

The  differential  equations  governing  the  growth  of  radical  concentrations  were  then  set  up 
in  Reference  2  for  the  induction  . period  and  a  particular  set  of  solutions  were  assumed,  of  the 
form 

Cj  =  A|  ext  i  *  OH,  H,  0,  H02 

For  the  case  of  interest  below. the  second  explosion  limit  where  2k2  >  k# 

^  >,^(21,^  t,CM)CB|Co2  _ 

S  Sch*  +  (VS)  k«  cm2c02cm+[  (hz+ k$ CM) Cq^]  X +Xz 
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For  the  case  where 

*6  CM  >  2kZ 

x  «  2k2  *11  ^2 

k6  CM  -  2k2 

The  end  of  the  induction  period  was  then  defined  to  be  when  the  concentration  of  OH  reached 
-6 

10-  mole/liter,  with  the  result 

t  =  <b  -  £n  P)/X  (is) 

The  value  of  the  constant  B  was  taken  to  be  25  to  correlate  high  temperature  Ignition  delay 
data  where  2k2>k,cM.  For  longer  delays  where  k,  cM>2kl<  a  vftlue  of  B  *  3.8  was  found  to  cor¬ 
relate  the  experimental  data.  Equation  15  was  then  used  to  compute  values  for  the  ignition 
delay  in  the  temperature  range  1500-2000®R  at  pressures  of  15  and  30  psia.  The  results  are 
shown  in  Figure  7. 

The  theoretical  analysis  of  molecular  reaction  rates  of  Reference  5,  which  has  been  used 
with  good  success  in  predicting  reaction  rates,  indicates  that  recombination  reactions,  such 

-0  4 

as  reaction  6.  have  rate  constants  proportional  to  T  *  .  This  temperature  dependence  has 
been  experimentally  undetectable  due  to  inaccuracies  In  the  experimental  data,  although  It  was 
believed  that  a  temperature  dependence  should  exist.  The  value  of  kg  was  then  computed  at 

300°K  from  the  equation  of  Reference  2.  Using  this  value  for  kg  at  3001^  and  assuming  the  rate 
to  be  proportional  to  T~0,4,  a  new  effective  collision  frequency  was  computed  giving 

k,  =  6.24  x  lO10  (XH  +  0.35  X„t  +  0.43  XNj  + - )  T-0,4 

With  this  new  value  for  kg.  the  ignition  delay  calculations  were  repeated.  The  results  are 
shown  in  Figure  8. 
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Figure  8.  Ignition  Delay  of  Hg-Alr  Versus  Temperature 
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SECTION  IV 

EXPERIMENTAL  EQUIPMENT  AND  PROCEDURE 


The  shock  tube  used  in  this  study  consisted  of  a  driver  and  a  driven  section,  each  constructed 
of  3-inch  stainless  steel  pipe,  20  feet  long,  with  no  internal  machining.  Attached  to  the  end  of 
the  driven  section  was  a  2-foot  test  section  containing  a  Kistler  Model  603  quart*,  pressure 
transducer  for  monitoring  the  reflected  shock  pressure  and  two  1/2-lnch-diameter  quartz 
windows  for  monitoring  the  OH  radiation  during  the  combustion  process.  The  quartz  windows 
and  pressure  transducer  were  located  4  inches  from  the  end  of  the  test  section.  At  a  distance 
of  1.905  feet  up&tream.  another  pressure  transducer  was  installed  to  trigger  the  electronic 
equipment.  The  time  required  for  the  shock  wave  to  traverse  the  distance  between  the  two 
pressure  transducers  was  recorded  on  two  Atec  electronic  counters. 

The  light  output  from  the  combustion  process,  initiated  by  the  reflected  shock,  was  passed 
through  an  Edmund’s  Scientific  Company  diffraction  grating  monochrometer  set  at  306  milli¬ 
microns.  The  light  intensity  from  the  monochrometer  was  then  measured  by  an  RCA  photo¬ 
multiplier  tube.  This  data,  together  with  the  pressure  measurements  from  the  transducer, 
were  recorded  on  three  Tektronic  535  oscilloscopes. 

The  pressure  transducers  were  mounted  in  the  tube  by  means  of  a  shock  mounting  system 
employing  neoprene  “O”  rings  to  completely  isolate  the  transducers  from  the  tube.  The  mount¬ 
ing  system  was  required  to  reduce  the  vibrations  picked  up  by  the  transducers  from  the  tube. 
The  signals  from  the  transducers  were  amplified  by  two  Kistler  566  electrostatic  charge 
amplifiers  and  recorded  on  the  Tektronic  oscilloscopes. 

The  double-diaphragm  technique  was  used  for  rupturing  the  diaphragms.  Two  Mylar  dia¬ 
phragms  of  0.001,  0.002,  or  0.003  inch  thickness,  depending  on  the  final  pressure  in  the  driver 
tube,  were  separated  by  a  1-inch  metal  plate.  The  driver  section  and  the  section  between  the 
diaphragms  were  loaded  to  one-half  the  final  pressure.  The  driver  was  then  pressurized  to 
its  predetermined  final  pressure,  and  the  center  section  between  the  two  diaphragms  was 
vented  to  a  vacuum  tank,  allowing  the  first  diaphragm  to  be  ruptured  and  then  the  second. 

The  hydrogen-air  mixture  was  premixed  in  stoichiometric  proportions  in  a  clean  gas  bottle 
and  filled  to  a  pressure  of  about  100  psla.  The  mixture  was  Introduced  into  the  driven  tube  to 
the  desired  pressure  after  the  tube  had  first  been  evacuated  to  about  200  microns.  This  pres¬ 
sure  was  measured  with  a  Wallace  and  Ticrnan  0-200  mm  Hg  absolute  pressure  gage. 

For  the  tailored  mode  of  operation,  the  nitrogen  pressure  in  the  driven  tube  was  measured 
with  a  0-200  Inches  of  mercury  Kollsman  gage.  The  final  driver  pressure  wa3  recorded  on  a 
0-100  psia  Brown  recorder. 
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SECTION  V 

DISCUSSION  OF  EXPERIMENTAL  DIFFICULTIES 


During  the  course  of  this  investigation,  many  difficulties  were  encountered  which  made  the 
evaluation  of  the  measurements  initially  impossible. 

The  first  of  these  difficulties  was  in  obtaining  proper  triggering  of  the  electronic  instru¬ 
mentation  as  the  shock  wave  passed,  which  was  necessary  to  make  accurate  shock  speed  mea¬ 
surements.  This  problem  was  at  first  attributed  to  high  humidity  since  the  room  in  which  the 
experiments  were  being  performed  was  not  climate -controlled.  Baking  the  transducers  and 
cables  did  seem  to  help,  but  only  temporarily.  The  oo-axial  cable  between  the  transducers  and 
charge-amplifiers  was  finally  removed  and  the  transducers  were  connected  to  the  charge 
amplifiers  with  short  lengths  of  Kistler  low-noise  cable.  This  change  completely  eliminated 
all  difficulties  in  triggering  and  shock  speed  measurement. 

The  next  problem  was  not  as  vital  to  the  completion  of  the  study  as  the  first  but  was  equally 
perplexing.  As  shown  in  Figure  9,  all  oscilloscope  traces  of  the  reflected  shock  pressure  at 
tailored  conditions  exhibited  a  hump  in  the  pressure  trace  approximately  2  milliseconds  after 
the  reflected  shock  passed.  Although  the  pressure  rise  was  only  on  the  order  of  10%,  we  needed 
to  know  whether  this  disturbance  sms  caused  by  the  specific  shock  t’^oe  being  used  in  this  in¬ 
vestigation  or  if  it  was  characteristic  of  all  shock  tubes.  No  mention  of  an  equivalent  problem 
was  found  in  the  literature.  A  great  deal  of  effort,  . therefore,  was  expended  in  an  attempt  to 
eliminate  this  problem,  but  with  no  success.  It  was  noticed  after  some  time  that  the  position 
of  this  hump  varied  noticeably  with  shock  velocity,  which,  of  necessity,  implies  that  the  dis¬ 
turbance  is  produced  from  the  contact  surface  as  foe  reflected  shock  wave  passes.  We  con¬ 
cluded  that  a  region  of  mixing  exists  at  foe  oontect  surface  between  the  driver  and  driven 
gases,  such  that,  no  matter  what  foe  shock  speed,  it  is  always  overtailored.  A  series  of  weak 
compression  waves  and,  subsequently,  expansion  waves  are  then  produced  as  foe  reflected 
shock  passes  through  this  region.  These  conclusions  were  later  confirmed  by  foe  results 
presented  in  Reference  6. 

Another  'problem  was  encountered  when  foe  ignition  delay  experiments  were  started.  The 
problem  encountered  at  this  point  proved  to  foe  the  most  difficult  of  foe  entire  experimental 
program.,  When  foe  hydrogen-air  mixture  was  introduced  into  the  driven  tube  and  the  re¬ 
flected  shock  was  used  to  bring  foe  temperature  and  pressure  of  foe  mixture  up  to  a  point 
where  spontaneous  combustion  could  occur,  foe  mixture  would  ignite  behind  foe  incident  shock, 
in  a  region  where  the  temperature  was  only  70p-800*R,  which  is  Well  below  the  explosion 
limits  for  hydrogen-air.  Observation  of  shock-heated  air  through  foe  quartz,  windows  showed 
foe  tube  to  contain  burning  particles.  The  tube  was  foen  completely  disassembled  and  cleaned, 
but  this  did  not  change  foe  behavior.  Filters  were  placed  on  all  gas  lines  entering  foe  tube,  but 
this  did  not  solve  the  problem.  Efforts  to  identify  foe  nature  of  foe  particles  were  unsuccessful. 
A  fine  dust  was  noticed  around  the  double  diaphragm  section  of  foe  shock  tube,  however,  and 
when  this  dust  was  introduced  around  foe  quartz  windows  of  foe  tube,  a  brilliant  flash  was 
observed  at  the  window  station  when  foe  tube  was  fired.  We  noticed  that  foe  arrangement  of 
the  double  diaphragm  section  was  such  that  small  metal  particles  from  foe  bars  holding  the 
double  diaphragm  were  scraped  off  and  entered  foe  tube  each  time  it  was  closed.  The  support 
bars  were  moved  to  the  side  of  foe  tube,  which  seemed  to  end  foe  problem  and  allowed  foe 
initial  ignition  delay  data  to  be  taken. 

Over  a  period  of  several  weeks,  ignition  delay  date  for  reflected  shock  pressures  of  15  and 
30  psia  and  temperatures  of  1580*R  to  1800®R  were  obtained.  A  sample  oscilloscope  trace  is 
shown  in  Figure  10  for  P  »  30  psia  and  T  «  1705*R.  The  results  of  these  experiments  are 
tabulated  in  Table  1  and  shown  in  Figures  11  and  12,  along  with  foe  analytical  predictions  for 


17 


> 


A  FAPL-TR-66-74 


the  two  values  of  the  rate  constant  kg.  Figure  13  shows  how  the  data  relates  to  the  correlation 

of  high-temperature  data  and  emphasizes  the  fact  that  the  high-temperature  correlation  cannot 
be  used  below  temperatures  of  1700  -  2000CR,  depending  upon  the  pressure  of  the  mixture. 


The  temperatures  for  the  data  shown  were  calculated  from  the  measured  shock  wave  velocity. 
Several  means  of  computing  the  temperatures  were  tried  to  see  If  the  scatter  in  the  data  could 
be  reduced,  such  as  computing  the  temperature  from  the  measured  reflected  shock  pressure, 
and  computing  a  temperature  time  history  of  the  ignition  period  from  the  pressure  trace  and 
time  averaging  the  curve  to  get  an  average  temperature.  Neither  of  these  methods  appeared 
to  reduce  the  scatter  of  the  data  significantly. 

The  data,  at  first  glance,  appear  to  be  quite  scattered,  but  all  but  8  of  the  43  experimental 

-0  4 

points  are  within  a  ±2%  temperature  band  of  the  analytical  prediction  for  kg  a  T  '  .  This 
would  seem  to  lend  support  to  the  adoption  of  the  reaction  rate 

k,  =  6.24  x  I010  (*Hi+  0.35  x^  +  0.43  XNi+ - )  T~a4 

This  was  a  very  limited  amount  of  data  on  which  to  base  such  a  conclusion,  however,  so 
additional  experiments  were  planned.  Since  there  is  less  slope  to  the  computed  ignition  delay 
curve  at  low  pressures,  the  obvious  approach  would  be  to  use  reflected  shock  pressures  of 
5  or  10  psia.  The  electr  jnlc  counters  required  a  signal  of  approximately  250  millivolts  for  re¬ 
liable  triggering  and  the  pressure  transducer  output  and  amplifier  gain  produced  only  39  mil¬ 
livolts  per  psi.  however,  so  it  was  impossible  to  proceed  with  these  lower  shock  pressures. 

The  next  series  of  experiments  was  then  planned  to  be  conducted  at  60  psia.  In  this  set  of 
experiments,  the  problem  of  ignition  well  below  the  ignition  limits  was  again  encountered, 
although  the  15  and  30  psia  data  were  still  repeatable.  Lowering  the  pressure  to  45  psia  pro¬ 
vided  no  improvement  nor  did  dismantling  and  cleaning  the  tube. 
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TABLE  I 

IGNITION  DELAY  OF  Hj-AIR 


(<f>  -D 


15 

psla 

30  psla 

T  (*R) 

ms 

T  (*R) 

ms 

1575 

3.25 

1677 

0.6 

1655 

0.3 

1677 

0.9 

1655 

0.3 

1702 

0.85 

1580 

0.595 

1705 

0.75 

1648 

0.51 

1796 

0.080 

1690 

0.41 

1759 

0.078 

1682 

0.41 

1735 

0.105 

1710 

0.225 

1725 

0.215 

1725 

0.27 

1740 

0.140 

1710 

0.25 

1688 

0.130 

1767 

0.11 

1700 

0.245 

1704 

0.355 

1720 

0.16 

1735 

0.205 

1680 

0.52 

1600 

0.7 

1692 

0.4 

1582 

1.85 

1655 

2.05 

1591 

1.6 

1645 

0.45 

1775 

0.12 

1655 

0.44 

1575 

2.4 

1665 

0.25 

1631 

0.48 

1665 

0.2 

1580 

1.25 

1640 

3.4 

1596 

1.5 

1652 

1.65 

1595 

1.35 

1635 

2.25 

EXPERIMENTAL 
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Figure  13,  Ignition  Delay  of  Stoichiometric  Hg-Air  Mixture  Versus  Temperature 
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SECTION  VI 
CONCLUSIONS 


A  shock  tube  was  constructed  and  ignition  delay  data  at  IS  and  30  psia  successfully  ootained 
near  the  second  explosion  limit  of  hydrogen  -  air.  Since  the  data  were  obtained  for  a  stoichio¬ 
metric  mixture  with  no  diluent,  the  data  do  not  require  extrapolation  for  use  in  air-breathing 
engine  calculations.  Also,  the  data  emphasise  the  fallacy  of  extrapolating  high  temperature 
ignition  delay  data  for  hydrogen  -  air  for  use  at  conditions  near  the  second  explosion  limit 
In  addition,  the  analytical  technique  of  Reference  2,  when  applied  with  a  modified  reaction 
rate  for  the  H02  recombination  reaction,  does  predict  the  experimental  results  well  within 

the  accuracy  of  the  data. 

Many  problems  were  brought  to  light  and  resolved,  which  should  make  further  studies  much 
less  difficult  and  time  consuming.  A  severe  problem  was  encountered  because  the  tube  was 
constructed  of  stainless  steel  pipe  of  unknown  origin,  and  there  were  many  crevices- in  which 
foreign  matter  oould  accumulate.  Future  tubes  will  be  constructed  of  new  tubing  which  has 
been  carefully  reamed  and  cleaned,  and  which  will  have  a  minimum  of  instrumentation  parts. 

The  results  of  this  study  suggest  that  additional  work  would  be  desirable  in  the  following 
areas:  (1)  data  obtained  at  higher  and  lower  pressure  than  the  present  study  to,. hopefully, 
provide  additional  support  for  the  recommended  value  of  kg  and  (2)  addition  of  various  diluents 

to  the  hydrogen  -  air  mixture  to  see  if  effect  on  foe  H02  reaction  is  as  given  in  Reference  2. 

One  other  area  to  be  investigated  is  the  effect  of  catalysts  on  the  ignition  delay.  It  is  known 
that  the  small  amounts  of  N02  reduce  the  ignition  delay  at  low  temperatures  by  at  least  an 

order  of  magnitude;  an  investigation  should  be  conducted  to  see  if  this  effect  is  due  to  the 
NO  released  by  the  dissociation  of  the  N02  and  reacting  with  the  Hg  or  merely  the  atomic 

oxygen  released  by  fee  dissociation  of  fee  N02.  In  addition,  fee  use  of  a  fine  dust  as  an  ignition 

source  at  very  low  temperatures  should  be  investigated. 
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